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ABSTRACT 



Aims. Supernova remnant (SNR) RX J0852. 0-4622 is one of the youngest and is most likely the closest among known 
galactic supernova remnants (SNRs). It was detected in X-rays, the ''''Ti 7-line, and radio. We obtain and analyze 
medium-resolution spectra of 14 stars in the direction towards the SNR RX J0852. 0-4622 in an attempt to detect broad 
absorption lines of unshocked ejecta against background stars. 

Methods. Spectral synthesis is performed for all the stars in the wavelength range of 3740-4020A to extract the broad 
absorption lines of Ca II related to the SNR RX J0852.0-4622. 

Results. We do not detect any broad absorption line and place a Scr upper limit on the relative depths of <0.04 for the 
broad Ca II absorption produced by the SNR. We detect narrow low and high velocity absorption components of Ca II. 
High velocity IVlsaI ~ 100 — 140km s~^ components are attributed to radiative shocks in clouds engulfed by the old 
Vela SNR. The upper limit to the absorption line strength combined with the width and flux of the ^^Ti 7-ray line 
1.16 MeV lead us to conclude that SNR RX J0852. 0-4622 was probably produced by an energetic SN Ic explosion. 



Key words. Line: formation - Stars: fundamental parameters - 
RX J0852. 0-4622 - ISM: supernova remnants 
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1. Introduction 

Young supernova remnants (SNRs) at the deceleration 
stage generally provide us with an opportunity to probe the 
supernova (SN) ejecta not yet polluted by the circumstel- 
lar matter (CSM). The ejecta composition and structure 
are usually studied by analyzing the X-ray spectra ema- 
nating from the reverse shock. In some young SNRs related 
to core-collapse SNe, e.g., Cas A, the optical emission of 
undecelerated ejecta clumps are observed when they pene- 
trate the post-shock layer, in which case they are powered 
by slow radiative shocks. Unshocked ejecta material is cold 
and does not radiate, thus remaining invisible in emission. 

However, in rare cases the unshocked ejecta of SNR 
can be observed in resonance absorption lines against 
a background light source seen through the SNR shell. 
Among Galactic SNR, this method has been applied suc- 
cessfully only to SN 1006, where the ejecta was observed in 
the ultraviolet a bsorption lines against the spectra of hot 
stars and QSOs (Hamilto n et al.lll997HWinkler et al.ll2005l: 
[Hamilton ct al. 2007). There is also one extragalactic SNR 
detected in absorption lines: the remnant of SN 1885 in 
M 31. It was first detected by groun d-based imaging in Fe I 
3860 A band against the M 31 bulge (iFesen et al.lll989l) and 
afterw ards by HST imaging and spectroscopy (^ Fesen et all 
Il999f ). The optical spectrum of this SNR contains strong 
Fe I, Ca II, and Ca I broad resonance absorption lines pro- 
duced by the unshocked ejecta. 



RX J0852.0-4622 (Vela Jr., G266.2-1.2) is a young galac- 
tic SNR detected by m eans of its emission in hard X- 
(Aschenbachl | 1998D . t he ^^Ti 1.16 MeV 7-ra y line 



rays 



(jlvudin et all I1998D. radio (iDuncan k GreenI \200(i ). and 



TeV ^-ravs (jAharonian et al.ll2005h . Vela Jr. with a diam- 
eter of ^ 2° is superimposed on the eastern part of the 
well-known old Vela SNR. The age and distance of Vela Jr. 
are estimated to 700 ± 150 yr and ~ 200 p c, respectively 
(jAschenbach et al.lll999l : iBamba et al.ll2005h . The age, dis- 
tance, and angular radius imply an average expansion ve- 
locity of 5000km s~^. 

Using XMM-Newton images, iKatsuda et al.l 



measured the proper motion of the bright NW rim of RX 
J0852. 0-4622. The derived value turns out to be about 
5 times lower than the predicted average expansion rate 
of Vela Jr. On the b asis of the measured proper motion, 
IKatsuda et al.l (|2008D estimate the age of Vela Jr. to be 
1700-4300 years and its distance to be ~750 pc. Although 
the conclusion could be hampered by a possible interaction 
with the recently encountered dense interstellar medium, 
the possibility of a large age cannot presently be fully ruled 
out. 

Given the uncertainty in the age issue, we should care- 
fully study the implications of a young age. In this respect, 
it is tempting to consider the unshocked ejecta of Vela Jr. 
using absorption spectroscopy against background stars. 
Our preliminary estimates showed that broad Ca II ab- 
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Table 1. List of the observed stars 



N 


Star 


RA 


Decl 


V 


Sp 


B — V) 


P 


S/N 


O" 






(2000.0) 


(mag) 




(mag) 








1 


HD75309 


08 47 28.0 


-46 27 04 


7.84 


B2Ib/II 


0.01 


0.79 


270 


0.008 


2 


HD75820 


08 50 26.0 


-46 14 53 


8.64 


B9V 


-0.02 


0.27 


260 


0.007 


3 


HD75873 


08 50 48.8 


-46 18 36 


8.10 


A3II/III 


0.38 


0.20 


160 


0.010 


4 


HD75955 


08 51 26.0 


-45 37 23 


7.73 


B9V 


-0.01 


0.67 


240 


0.005 


5 


HD75968 


08 51 32.8 


-46 36 36 


8.14 


B9III/IV 


-0.12 


0.33 


200 


0.006 


6 


HD76060 


08 52 02.4 


-46 17 20 


7.88 


B8IV/V 


-0.09 


0.02 


350 


0.005 


1 


tliJ < D0o9 


08 55 23.0 


-46 53 28 


Q Q /I 




-U.Uo 


U.o4 


Z i U 


n Ann 

u.uuy 


8 


HD76649 


08 55 50.4 


-46 20 30 


8.33 


B7II/III 


0.14 


0.67 


150 


0.008 


9 


HD76744 


08 56 18.2 


-46 19 57 


8.69 


AOV 


0.08 


0.75 


200 


0.009 


10 


CD-454590 


08 49 35.5 


-46 23 18 


9.58 


B5 


0.20 


0.42 


140 


0.010 


11 


CD-454606 


08 50 15.0 


-45 31 22 


8.96 


B0.5V 


0.38 


0.82 


180 


0.009 


12 


CD-454645 


08 51 34.9 


-46 09 54 


10.32 


AO 


0.20 


0.14 


90 


0.012 


13 


CD-454676 


08 53 22.0 


-46 02 09 


8.93 


B0.5III 


0.77 


0.35 


140 


0.014 


14 


CD-464666 


08 50 44.3 


-46 38 11 


9.81 


AOII 


0.60 


0.40 


90 


0.009 



sorption lines might be observed. Here we report results of 
spectral observations and analysis of the spectra of distant 
stars across the Vela Jr. In Sect. 2, we describe the ob- 
servations and data reduction. The results of the spectral 
synthesis and extraction of broad and interstellar Ca II ab- 
sorption are presented in Sect. 3. We fail identify any broad 
absorption lines and the implications of this are discussed 
in Sect. 4. 

2. Observations and data reduction 

Spectra were obtained on the ESO 3.6-m telescope 
NTT (program 080.D-0012(A) PI: A.F.Iyudin) using the 
echelle spectrograph EMMI in BLMD mode of medium- 
dispersion spectroscopy with holographic grating =1^11 
(3000 grooves/mm and maximum light reflectivity at 
Ab;aze=3500A). The disper sion is O.lSA/pix, the resolving 
power is A/AA=9000 for a slit width of 1.02" and a regis- 
tered wavelength band of ~1500A. The signal-to-noise ratio 
(S/N) is between 90 and 350. For each star, two overlapped 
spectra in the range of 3740-4021A were obtained with cen- 



ters at 3820A and 3945A. This spectral region encompasses 
resonance spectral lines of Fe I (3860A) and Ca II (3933A, 
3968 A). 

Among the observed stars, ten are of B-type and four of 
A- type (Table [J) with magnitudes 7™ < V < 11™. Tabled] 
contains the star number, star HD/CD name, right ascen- 
sion, declination, F-magnitude, spectral type, B — V color 
index, relative impact parameter p, and signal-to-noise ra- 
tio {S/N), and Icr, which represents the root mean square 
error in the spectral fit (see Sect. 3). The impact parame- 
ter is defined as p = 9 / p, where 9 is the angular distance 
of the star from the SNR center and p = 1° is the SNR 
angular radius. The SNR center presumably coin cides with 
the n eutron star candidate AX J0851. 9-4617 f|Slane et al.l 
I2"00l] ). which has coordinates (RA,Dec)2ooo = (08'' 51"" 57^ 
—46° 17'. 4). The preliminary estimated distances of selected 
stars are in the range between 240-2000 pc, so at least some 
of the stars lie behind Vela Jr. Positions of all stars of the 
program across the Vela Jr. X-ray image (Rosat All Sky 
Survey data with energies ab ove 1.3 keV) and TeV 7-ray 
image (EE aroman et al.ll2005ll are shown in Fig. [1] 




09:00 



08:50 



Fig. 1. Positions of the stars across the SNR images taken 
in TeV 7-rays (grayscale) and in hard X-rays (contours). 
The star numbers correspond to those in the first column 
of Table [TJ The cross is the position of the neutron star 
AX J0851. 9-4617. The size of symbols qualitatively reflects 
the distance with the largest symbol being closest star. 



Preliminary processing of the spectral CCD images, 
extracting the spectra, and performing wavelength cal- 
ibration using ThAr-lamp data were all completed by 
employing the EMMI native emmi-quickred script of 
MIDAS; the atmosphere and interstellar extinction was 
also taken into account for each individual star. The 
spectra are flux calibrated with the help of spectra of 
the standard star HD 60753 taken in the same set 
of observations. This star has a calibrated flux in the 
UV band only. The v isual flux is synthesized applying 
Kuru cz ATLAS9 code (iKuruczl [19931) and SynthVb pack- 
age (|Tsvmbal et al.l I2003D . AT stellar atmosphere model 
is calculated for solar metallicity, effective temperature 
Te//=16200 K ((7=150 K), gravity log 5=3.56 (cr=0.05), 
rotation velocity y sini=24 km s~^{a=2 km s~^), and in- 
terstellar absorption Ay = 0.25 (cr=0.03). The derived re- 
sponse function is used to produce flux-calibrated stellar 
spectra. One should emphasize that this procedure provides 
us with the precise relative flux across a wide wavelength 
range but not the absolute flux; the latter, in fact, is not 
needed for our purposes. 
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Table 2. Parameters of the stars and distances calculated by spectral method and from Hipparcos parallax 



3 



Star 


Teff 




log g 


V sin i 


Vr 


{B 


- V)o 


Av 


M 


dsp 


d-Hip 




(K) 






(km s~^) 


(km s"') 


(mag) 


(mag) 


{Mq) 


(pc) 


(pc) 


HD75309 


26500±400 


3 


60±0.10 


210 


43 




-0.25 


0.8 


16 


1900±300 




HD75820 


11400±200 


4 


OOiO.lO 


200 


31 




-0.10 


0.2 


3.2 


470±100 


505±184 


HD75873 


8900±200 


2 


50±0.05 


15 


41 




0.01 


1.2 


6 


1400±200 




HD75955 


10400±200 


3 


85±0.10 


190 


24 




-0.07 


0.2 


3.0 


320± 70 


262± 35 


HD75968 


12250±150 


3 


86±0.05 


80 


35 




-0.11 


0.0 


3.8 


570±140 


719±254 


HD76060 


13400±200 


4 


lOiO.lO 


240 


36 




-0.13 


0.1 


3.6 


390± 90 


335± 63 


HD76589 


11800±200 


4 


lOiO.lO 


95 







-0.10 


0.1 


3.2 


390± 90 


240± 76 


HD76649 


13300±150 


3 


65±0.05 


33 


33 




-0.13 


0.8 


4.5 


640±110 




HD76744 


10500±200 


4 


20±0.10 


150 


-7 




-0.08 


0.5 


2.4 


270± 50 




CD-454590 


22400±400 


3 


60±0.10 


140 


31 




-0.22 


1.3 


11 


2400±300 




CD-454606 


29500±500 


3 


80±0.10 


240 


33 




-0.27 


2.0 


17 


1670±160 




CD-454645 


8400±200 


4 


35±0.10 


80 


35 




0.08 


0.4 


1.8 


330± 70 




CD-454676 


29000±500 


3 


70±0.10 


125 


20 




-0.27 


3.2 


18 


1080±150 




CD-464666 


10500±200 


2 


05±0.05 


30 


49 




-0.08 


2.1 


12 


5700±500 





3. Analysis of spectra 

The obtained stellar spectra were analyzed by applying 
a synthetic flux calculation employing the SynthVb and 
ATLAS9 codes to derive the stellar parameters Te/f, log g, 
and V sin i (Table [2]) . To determine the rotation velocity, 
we use a standard method based on Fourier transformation 
of profiles of weak spectral lines (|Carrolllll933l) . We adopted 
the limb darkening parameter e = 0.6. 

The parameters of stellar atmosphere models were 
found by fitting the synthesized profiles to six Balmer lines 
(H7-H12). The fitting procedure includes multiplication of 
the observed spectra by a fitting factor C(A), which is de- 
fined as the linear approximation to the ratio of synthetic 
to observed spectrum F^^^ /F°^^ . In this case, we do not 
introduce any nonlinear distortion into the observed spec- 
trum. The best fit is attained by minimizing the relative 
residual fiux between the observed and synthetic spectrum 



C(A)*F/ 



'obs 



psyn 



'Syn 



(1) 



where F^^^ is the continuum fiux of the synthetic spectrum. 
Using the relative flux permits us to avoid problems related 




K (A) 

Fig. 2. Observed spectrum of the star HD 75968, synthetic 
spectrum, and the residual spectrum with 3fT error box 



-a" 



400 



600 

,(pc) 



Fig. 3. Distance derived by the method of spectral parallax 
versus distance according Hipparcos parallax. 



to a poorly deflned observed continuum. An example of the 
flt quality is demonstrated by Fig.[31 The la statistical error 
in the spectral fit for individual stars are given in the last 
column of Table [TJ Any interstellar lines with an intensity 
exceeding the 3cr error, which lie in the range of 1.5-4.2%, 
should be visible in the residual spectra. As an example, 
the residual spectrum of star HD75968 (Fig. [5]) appeanrs to 
contain the Ca II interstellar absorption doublet. We believe 
that the systematic errors in producing residual spectra are 
smaller than the quoted 3cr statistical error. 

3.1. Interstellar reddening and distance estimate 

Interstellar reddening in the observed sp ectra is taken into 
account using extinction data reported bv lMathi^ (jiggoT ) . In 
the first step, we determine stellar parameters and calculate 
normal color indexes {B — V)q (fCastclli & Kurucz 2003|). 
We then use the observed {B — V) to derive color excess 
E{B - V) and V-band absorption, Ay = 3.1E{B -V). At 
the second step, we improve our determination of the stellar 
atmosphere parameters and recalculate the reddening. 

To determine distances, we employ a modified method 
of spectral parallaxes in which the stellar luminosity is de- 
rived from stellar evolutionary tracks as follows. Using stel- 
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lar parameters (Tefj and log g) and evolutionary tracks 
(|Schaller et al.lll992ll . we estimate the stellar mass and thus 
derive the bolometric luminosity. The absolute magnitude 
My is then determined using the bolometric correction 
taken from (|Bessell et al.l [T998() . The luminosity and dis- 
tance are thus determined from the standard formulae 

log L = -10.607 + log (M/Mq) + 41og Te// - log g , (2) 

Mv = 4.69 - 2.5 log L - BCy , (3) 

log dsp = 0.2(My - mv -5 + Av). (4) 

The normal color index {B — V)o, interstellar absorption 
Ay, stellar masses M, and distances are listed in Tabled) 
In the last column, we also indicate th e ava ilable distances 
according to HIPPARCOS parallaxes (Ivan Lccuwcn 2007). 
The distances determined by both methods agree within 
errors (Fig. [3]), which supports the reliability of distances 
obtained by the method of spectral parallax. 

3.2. Interstellar lines 

The relative residual spectra for all the stars are shown 
in Fig. m The Scr levels are also indicated for each star. 
Large values of 3cr are seen for late B and A type stars 
with low rotation speeds due to problems in fitting narrow 
stellar spectral lines using the accepted solar abundance 
of chemical elements. The spectra do not uncover broad 
absorption resonance lines of Ca II 3933A, 3968A, or Fe I 
3860A. In particular, the relative depth of the broad Ca II 
absorption (if any) produced by Vela Jr. is smaller than 
0.04 at the level of 3cr. We note, that the weak absorption 
at 3819A and 4009A in the hottest stars of our sample are 
related to helium lines, which are generally affected by non- 
LTE excitation and cannot be modeled reliably in the LTE 
approximation. 

With the exception of two stars (HD 75873 and CD- 
454645), all the spectra contain narrow unresolved interstel- 
lar Ca II lines; their heliocentric radial velocities are given 
in Table [31 The interstellar absorption lines can be divided 
into two major groups: low velocity |T^|^50 km s~^, and 
high velocity [V'l^lOO km s~^. Most stars have one compo- 
nent with a positive radial velocity of ^22-48 km s^^. The 
heliocentric velocity can be translated into a LSR velocity 
in this direction using the relation Vlsr = Vhei — 13 km s~^ . 
This suggests that the dominant population of interstellar 
clouds in this direction at the distances not exceeding 2 
kpc, are characterized by the positive LSR velocities '^9- 
35 km s^^. Two stars have negative low velocity compo- 
nents, of -13 km s~^ and -46 km s~^. 

Three stars have high velocity components: HD 75309 
(-fl53 km s-\ -92 km s^^), HD 76060 (-92 km s'^), 
and CD-454676 (-150 km s"^). These velocities are typi- 
cal of high- velocity interstellar Ca II absorption found ear - 
her in the direction of Vela SNR (jCha fc SembachI |2000D . 
Interestingly, the spectrum of the star HD 75309 from 
our program was also studied by Ch a fc SembachI ()2000f) . 
Benefitting from the spectrum's high resolution, these au- 
thors were able to find eight components including two high 
velocity components, -1-136 km s~^ and -107 km s~^, and 
a strong low velocity component -1-20 km s~^. The corre- 
sponding interstellar absorption lines in Table [3] are shifted 
redward by « -1-16 km s~^, which reflects the systematic 
difference in radial velocity between the two sets of data. 



Table 3. Velocities of components of Ca II interstellar ab- 
sorption. 



Star 3933^ 3968^1 
A(A) V(kms~^) A(A) V(kms~^) 



J-J-L-' t OOkJU 


3Q35 67 


153 


3970 53 


156 




3934 1 3 


36 


3968 91 


33 




3933 06 


-46 


3967 84 
ouyj t .04: 


-48 




3932 46 


-92 


3967 23 


-94 




3933 98 


25 


3968 76 


22 


Hn7,S873 


3933 98 


25 


3968 76 


22 


J-J-J--' 1 OiJOO 


3933 9^1 


22 


3968 73 


19 


J-J-J--* 1 OiJ\J(j 


3933 99 


25 


3968 76 


22 


j-ij--' i \jyj\jyj 


OiJO'-ii. J- 4: 


■Hfi 
ou 


oyuo. ( u 


22 




3932.46 


-92 


3967.23 


-94 


HD76589 


3934.10 


33 


3968.88 


31 


HD76649 


3934.14 


36 


3968.76 


22 


HD76744 


3933.95 


22 


3968.73 


19 


CD-454590 


3934.10 


33 


3968.88 


31 


CD-454606 


3933.11 


34 


3968.88 


31 




3933.50 


-13 


3968.27 


-15 


CD-454645 


3934.29 


48 


3969.07 


45 


CD-454676 


3933.99 


25 


3968.76 


22 




3931.70 


-150 


3966.47 


-151 


CD-464666 


3934.20 


41 


3968.87 


30 



This may be partially related to the low resolution of our 
spectra. We studied other sources of errors but were unable 
to explain this disparity. 

At least one star, CD-454676, exhibits conspicuous 
CN absorption of electronic transitions R{0), -R(l), and 
P(l) with the wavelengths of 3873.994A, 3874.602A, and 
3875. 759A, respectively. The heliocentric radial velocity of 
these lines is -1-23 km s~^ {Vlsr = +10 km s^^), which 
is consistent with the radial velocity of Ca II interstellar 
lines in the same star (Table [21) • The equivalent width of 
-R(O) and i?(l) lines {W{0) = O.OlA and W{1) = O.O2A) 
can be used to estimate the excitation temperature of the 
rotational level J = 1 and the column density of CN re- 
siding on rotational levels J — and J = 1 in the weak 
line limit. Using available oscillator strengths of these tran- 
sitions (Gredelet aL ,200^), we obtain T — 6.8 K and 
column densities N{0) = 1.34 x 10" cm'^ and N{1) = 
1.16 X 10^'^ cm~^. Assuming a Boltzmann population of the 
J = 2 rotational level 7V(2) = 57V(0) exp(-16.325/r), we 
obtain the total CN column density -/V = 3.11 x 10^^ cm~^. 
These values are comparable to those of CN absor bers to- 
wards the Vela OB association (jGredel et al.ll2002l ). 

Two stars, HD 75873 and CD-454645, do not exhibit 
interstellar Ca II lines. In HD 75873 for which dgp— 1^00 pc 
and Ay— 1-2, the expected contribution of the interstellar 
line to the equivalent width should be about 20%. The ex- 
planation of the apparent absence of absorption is the low 
rotation velocity V smi=15 km s~^ and the high strength 
of stellar Ca II absorption. Both factors prevent us from dis- 
tinguishing interstellar lines in this case. The second star, 
CD-454645, at a distance of 330 pc is not expected to have 
strong interstellar Ca II lines. Given its very strong stellar 
Ca II absorption, the extraction of weak interstellar absorp- 
tion in this case is precluded. 
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Fig. 4. The relative residual spectra with dashed lines showing the 3cr statistical error boxes varying from ±0.015 to 
±0.04 from star to star (for Icr errors in each case, see Table [T]). The distance to each star is also indicated. 



Table 4. Adopted parameters of supernovae 



Parameters 


SN la 


SN IIP 


SN The 


SN Ic(h) 


M (Mq) 


1.4 


10 


3 


4 


E (lO^^erg) 


1.4 


1 


1.5 


20 



4. Discussion 

4.1. Narrow interstellar lines 

Most low velocity interstellar Ca II absorbers with Vlsr < 
20 km s~^ are most likely produced by local background 
clouds. This is also true for CN absorbers, which in terms 
of velocity coincide with the Ca II absorbers. The large scat- 
ter in velocities of between -50 and ±30 km s~^ exceeding 
the usual dispersion in cloud velocities of ~ 10 km s~^ sug- 
gests that at least some of these absorbers are related to 
clouds accelerated by either shock waves driven by wind 
bubbles of hot stars or old SNR. This conjecture is in ac- 
cord with results of observatio ns of interstellar ultrav iolet 
O I and Si II absorption lines (jWallerstein et al.lll995D to- 
wards the Vela SNR. Some of these lines arise from ex- 
cited fine-structure levels in the Vela SNR direction and 
indicate the high pressu re of clouds, p ^ 10^^° dyn cm~^ 
(jWallerstein et al.lll995[ ). which is two orders of magnitude 
higher than the average pressure in the interstellar medium 
(ISM). 

The high velocity clouds with |Vi,sj^ | ^ 100—150 km s~^ 
are probably related to the interstellar clouds shocked by 



the expanding Vela SNR. Similar high velocity gas was 
observed in Ca II lines and ultraviolet lines correspond- 
i ng to different ions, iiicludin g C I, O I, Mg I, and Mg II 
([Jenkins fc WallersteinI [19951 ) . and attributed to radiative 
shocks driven by the Vela SNR into interstellar clouds of 
density ~ 10 cm~^. Interestingly, the closest star with a 
high velocity interstellar component, HD 76060, lies at the 
distance 335±63 pc. This immediately provides an upper 
limit to the distance of the Vela SNR of ~335±63 pc, 
which is consistent w ith the Vela pulsar distance of 294 pc 
([Caraveo et al.ll200H) . 



4.2. Broad absorption related to Vela Jr. 

The absence of Fe I and Ca II broad absorption lines in stel- 
lar spectra towards Vela Jr. requires explanation. As we will 
see below, singly ionized metals should be more abundantly 
present in the Vela Jr. ejecta, and therefore absorption by 
neutral iron, with its relatively low ionization fraction and 
low value of the oscillator strength, should be significantly 
weaker than Ca II absorption. We therefore concentrate on 
the absence of broad Ca II lines. At least four possibilities 
are conceivable: Vela Jr. is farther than the most distant 
star in our sample; the SNR is much older; and the Ca II 
ionization fraction is small, i.e., Ca resides predominantly 
in the Ca I or in the Ca III ionization state. Discussion of 
these possibilities requires modeling the broad Ca II ab- 
sorptions expected at the given age for the remnants of the 
different SN types. 



A.F.Iyudin et al.: Broad absorption lines in the field of SNR Vela Jr. 




T3 



~1 ' \ 

CD-464666 



CD-454676 



HD76649 



HD75968 



HD75820 



HD76589 



HD76060 



CD-454645 



HD75955 



HD76744 



-400 



-200 





V(km/s) 



200 



400 



Fig. 5. The relative residual spectra for the stars. The dis- 
tance increases upward. 

4.2.1. Broad Ca II absorption for different supernova types 

The unshocked ejecta expands freely, i.e., the expansion 
law at a given age t is v — r/t. To describe this we will use 
cylindrical coordinates (z, p, 0) with z-axis coincident with 
the line of sight directed towards the center of SNR. The 
absorption produced by the scattering of the background 
stellar radiation in Call 3933, 3968 A lines at the radial ve- 
locity Vz is determined by a Sobolev optical depth in the 
resonant plane z — v^t along the line of sight of impact 
parameter p (assuming that the ejecta is spherically sym- 
metric) 



t{z,p,( 



0.0265/i2Ai2nit, 



(5) 



where the multipliers in the right-hand side in order are 
oscillator strength, wavelength, Ca II number density at the 
given radius r — (z^ + p'^Y^^, and the SNR age t. 

For a given concentration of Ca determined by the den- 
sity and Ca abundance, the line strength depends on the 
ionization fraction of Call. At the early phase of the ejecta 
expansion, at the time of t ~ 2 yr after SN explosion, the 
calcium ionization in ejecta of any type SN is controlled 
primarily by the ionization loss of fast electrons (Compton 
electrons and positrons) produced by the radioactive decay 
chain ^^Ni - ^^Co - ^^Fe and radiative recombination. At 
the stage of t > 3 yr, spectra of SNe of different types are 
dominated by the emission lines of singly ionized metals, 
which indicates that singly ionized metals dominate. This 
is supported by numerical models of ionization and thermal 
balance of ejecta powered by the radio active decay of ^^Co 
for SN la (|Axelrodlll98"ol ) and SN IIP (|Kozma fc Franssonl 
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Fig. 6. Absorption profile of Call doublet expected in the 
stellar spectrum of different progenitors of Vela Jr. Shown 
are cases of impact parameter equal to 0, 0.5, and 0.8. The 
strongest absorption always corresponds to zero impact pa- 
rameter. 



Il998[ ). At the later stages t ^ 10 yr, the ionization is domi- 
nated by positrons from ^^Ti decay. Our estimate indicates 
that even a maximal expected mass lO"** Mq of ^''Ti is 
insufficient to maintain the high ionization of Ca. At later 
stages, therefore, recombination dominates and Ca may be- 
come mostly neutral. 

However, at the SNR age of i ~ 10^ yr the characteristic 
recombination time of Ca is larger than the expansion time. 
At this stage, the ioni zation by the starlight may become es- 
sential. For example. iFesen et al.l (|1999 ) demonstrate that 
calcium in the ejecta of SN 1885 in M31 can be efficiently 
ionized by the bulge starlight within the ionization time of 
about 10 yr. In the case of Vela Jr. we use the model of the 
starlight spec trum in the Galact ic plane at the radius of 7.5 
kpc given bv lPorter et al.l (1200 6). and the photoionization 
cross-sections of lVerner et al.l (1199 61. The found photoion- 
ization time is ~ 180 yr for Cal and ~ 3 x 10^ yr for Call. 
At the age of ^ 700 yr, we thus expect that Ca in Vela Jr. 
should be singly ionized. 

The ionization of Ca II by X-rays from the reverse shock, 
and by accelerated protons may also play a role. The X-rays 
ionize metals from K and L shells, and photoelectrons then 



ionize Call. For the observed X-ray flux /x 



10- 



-10 



erg 



cm~ s~ in the 0.5-10 keV band, the characteristic pho- 
toionization time for Ca II at the SNR age of 700 yr is found 
to be ~ 10^ yr for SN la, and this process is thus negligible. 
Ionization by relativistic protons accelerated in the shock 
wave can be estimated by assuming an average efficiency of 
cosmic ray acceleration per SN of 10% and the shock wave 
energy of ^ 10^^ erg. We find then that the ionization time 
for Call at the age of 700 yr is ~ 2 x 10'^ yr, which is larger 
than the SNR age. We thus conclude that cosmic rays for 
the adopted acceleration efficiency essentially cannot ionize 
Call, so all the calcium in the unshocked ejecta of Vela Jr. 
is expected to remain in Call. 

The predicted profiles of the Ca II doublet at the age of 
700 yr for different varieties of SNe are shown in Fig. [6] as- 
suming that all the calcium is in the Ca II state. We assume 
that in SN IIP and SN Ibc the Ca abundance is solar, while 
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Fig. 7. Age-distance relations provided by ''^Ti mass 
(thick solid lines) and radius of the supernova remnant. 
The radius is calculated for SN IIP (dotted line), SN la 
(thin solid line), SN Ibc (short-dashed line), SN Ic(h) with 
35 Mq progenitor (long-dash line), and 60 Mq progenitor 
(dashed-dott e d lin e). Cross shows parameters derived by 
llvudin et all (|l998D . 

for SN la we assume that the Ca/Fe ratio by mass is solar, 
while the total mass of iron in the ejecta is 0.6 Mq. Ejecta 
parameters for different SNe are given in Table |31 Apart 
from SN la, SN IIP, and SN Ibc, we also consider energetic 
SN Ic, so-called hypernovae, which is designated hereafter 
as SN Ic(h). The boundary velocity of the unshocked ejecta 
is taken to be 5000 km s~^ in accordance with the distance 
of 200 pc and the age of 700 yr. The density distributions 
p(v) in the unshocked ejecta are assumed to be exponential 
for compact pre-SNe and to form a plateau with the outer 
power law p oc v~'^ for SN IIP. The plotted profiles are 
computed for three values of impact parameter in units of 
the angular radius: 0, 0.5, and 0.8. The absorption is pre- 
dicted to be deep for all the impact parameters in the case 
of SN la, rather deep for SN IIP, very weak for SN Ibc, and 
negligible (relative depth < 0.006) in the case of SN Ic(h). 
If the age and distance of Vela Jr. are close to the values 
adopted above, the progenitor would be unlikely of type 
SN la or SN IIP; instead an association of the SNR with a 
SN Ibc or SN Ic(h) is quite plausible. 

4.2.2. How far away might Vela Jr. be? 

We now relax arguments used earlier to con s train the age 
and distance of Vela Jr. (jAschenbachI 119981 : llvudin et al.l 
I1998D and check whether Vela Jr. lies at a very large dis- 
tance, > 1 kpc, beyond any star in our sample. 

For a given 44^j mass, a combination of age and distance 
is constrained by the observed flux in the gamma-ray line 
1.16 MeV. Nucleosynthesis mode ls predict a productio n of 
(1-5) X 10-5 Mq of "^^Ti in SN la (llwamoto et al.lll998ll and 
10-5 - 10-'' Mq in core-collapse SNe (jWooslev fc Weave^ 
Il995f) . An independent estimate of the '''*Ti yield per SN 
can be obtained by assuming that almost all the ^"^Ca 
is produced as '^''Ti bo th in SN la (llwamoto et al.l [l999l ) 
and core-collapse SNe (jWooslev fc Weaver! I1995D . The so- 
lar mass ratio of ''^Ca to ^^Fe is 10-^, which means that 
SN IIP and SN Ibc producing 0.05 - 0.1 of ^^Ni per SN 
should eject 5 x 10-^ - lO""* Mq of ''"'Ti, while SN Ic(h) 
producing, in a similar way to SN 1998bw, up to 0.5 M© 



of 56 Ni (jlwamoto et al.l Il998l ) is able to eject as much as 
5 X 10-* Mq of ''^Ti. We therefore, expect, that the mass 
of "''^Ti ejected by SNe of different types lies in the range 
5 X 10-5 — 5 X 10-* Mq. The corresponding relationship 
between the age and distance suggested by the observed 
flux of the 1.16 MeV line 3.8 x 10-5 cm^^ s^^ is shown in 
Fig. [Tjfor the two extreme values of ejected **Ti mass. 

Deceleration of supernova ejecta in the interstellar 
medium provides us with another relation between the dis- 
tance and age for a given choice of ejecta parameters, ISM 
density, and angular radius of Vela Jr. We compute the 
interacti on of ejecta wit h the ISM in the thin shell approx- 
imation (lChevaliedll982i) assuming typical ejecta mass and 
energy (Tabled]). In the case of SN IIP, the adopted hy- 
drogen number density of ISM is 0.3 cm-^, which is the 
average density of the warm neutral medium (WNM). The 
latter comprises about 80% of the ISM mass ([Wolfire et al.l 
1995). For SN la apart from WNM, we also consider the 
ISM in the form of a hot ionized medium (HIM) of density 
0.003 cm-'^. This ga s occupies about 50-60% of the volume 
(jWolfire et al.llT995l) . As in the case of SN Ibc and SN Ic(h), 
they explode in the ISM modified by the fast main-sequence 
wind, slow red supergiant wind, and the Wolf-Rayet wind. 
We consider 35 Mq and 60 Mq as template progenitor 
stars; both cases were explored by iGarcia-Segura et all 
(1996). According to these results the pre-SN in the 35 Mq 
case is imbedded in a hot bubble of the uniform density of 
^ 0.003 cm-"^ with a radius of 18 pc surrounded by a dense 
cool shell of total mass ^ 20 Mq. For the 60 Mq progenitor, 
the bubble density is ~ 0.001 cm-^ and its radius is 50 pc. 
We note in passing that a model of Vela Jr. (RXJ0852.0- 
4622) taken to be of the SNII/ Ib type exploded in a wind 
blown cavity was considered by iBerezhko et al.l ()2009l) . 

The age-distance relations for all the discussed cases 
are shown in Fig. jT] The SN la exploded in the HIM phase 
shows almost the same age-distance relation as SN Ibc and 
is therefore not shown in this figure. This diagn ostic plot 
is similar to that used bv lChen fc Gehrelsl ()1999t ). The es- 
sential difference, however, is that they used a set of arbi- 
trary expansion velocities of the swept-up shell, while we 
calculate the evolution of the shell radius for each type of 
SN. For a given age, the minimal distance corresponds to a 
SN IIP expanding in the WNM phase, while the maximal 
distance corresponds to a SN Ic(h) with a 60 Mq progen- 
itor. In combination with the curves, these two cases 
imply the allowed ranges of 450-900 yr and 150-1000 pc for 
the age and distance of Vela Jr., respectively. The major 
result of this plot is that the distance of Vela Jr. cannot ex- 
ceed 1 kpc. We thus conclude that at least several stars in 
our sample (Tabled]) are behind the SNR. This permits us 
to disregard the explanation of the absence of broad Ca II 
absorption being the result of the large distance to Vela Jr. 

4.2.3. Was the progenitor of Vela Jr. a hypernova? 

The absence of broad Call absorption in the spectra of 
stars at distances > 1 kpc suggests that the SNR progeni- 
tor was either of SN Ibc or SN Ic(h) because only for these 
SNe is the expected absorption weak and possibly unde- 
tected (Fig.O. To distinguish between these two SN possi- 
bilities, one should take into account the intrinsic width of 
the 1.16 MeV line. 

In the case of SN Ibc at the age of 650 yr, the expected 
profile of the **Ti line (cf. Fig. [7]) convolved with the instru- 
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mental profile (cr = 45 keV) is found to be too narrow com- 
pared with the observed one (Fig. |8^), even assuming ho- 
mogeneous mixing of 44q^j velocity of 10 000 km s~^. 
The observed broad profile of the 1.16 MeV line implies that 
most of the ejecta mass consisting of ''''^Ti has significantly 
larger velocities. However, the SN Ic(h) case with maximal 
expansion velocities of 31 300 km s~^at the age of ^ 500 yr 
(the case of 60 Mq) and spherically-symmetric distribution 
of ^^Ti homogeneously mixed to 31 000 km s~^ does not 
help resolve the ambiguity either (Fig. |S^). 

The solution to the line width problem might be found 
by taking into account that iron-peak elements are ejected 
by SN Ic(h) in the form of high velocity bipolar jets and as- 
suming that "^^Ti resides only in the outer parts of the jets. 
^^Ni-rich bipolar jets are predicted by the collapsar model 
(jMacFadven fc Woos lcv 1999) proposed for the hypernova 
SN 1998bw, and the jet-like structure is consistent w ith the 
spectr al line profiles (.Maeda et al..,2006,) . ,Maeda fc Nomotol 
(|2003) also predict an external location of '^'*Ti. In the 
1.16 MeV profile simulations, we assume that ^^Ti is ho- 
mogeneously distributed along the radius in the velocity 
range of 20000 — 31000 km s~^ within jets of opening angle 
60° and inclination angle 9. We took into account the light 
travel-time delay that produces the profile skewed towards 
red. Two cases are shown (Fig. |8}d) for the angle between 
the jet axis and the line of sight, 9 = 30° and 9 = 60°, 
both of which fit the data more closely than the spherically 
symmetric model. We therefore conclude that the hyper- 
nova model with the outer location of ^''^Ti in bipolar jets 
of SN Ic(h) is consistent both with the absence of the broad 
Ca II absorption, and the broad 1.16 MeV profile. 

A problem with the SN Ic(h) scenario is that the high 
ejecta velocity implied by this model infers a low ambient 
density, which seems to disagree with the baryonic origin of 
TeV gamma-ray emission from Vela Jr. Gamma ray produc- 
tion via pp-coUisions seems to be the pref erred model com- 
pared to the inverse Compto n mechanism (jAharonian et aP 
I2005t iBerezhko etal\ I20090 . An alternative may be pro- 
vided by assuming that we see the early stage of the in- 
teraction of the SNR with a dense environment that has 
not yet been affected by the previous expansion dynamics. 
This conjecture is in line with the low expansio n velocity 
found for the NW rim bv iKatsuda et all (|2008D . Another 
discomfort is related to the hypernova being a rare vari- 
ety of SNe that coniprises only about 1% of all SNe Ibc 
(jPodsiadlowski et"al]l2004[ ). Only high signal-to- noise spec- 
tral imaging of Vela Jr. in the 1.16 MeV line band with en- 
ergy resolution of < 40 keV and angular resolution of < 1° 
may confirm (or reject) the high velocities of ■*'*Ti and de- 
tect any jet-like (if any) structure of the distribution. 

Given the difficulties arising in the interpretation of 
data on Vela Jr., we should not ex clude out completel y 
the possibility that this SNR is older (|Katsuda et al.ll2008[) . 
However, only additional observations at different wave- 
length bands will be able to help pinpoint the age and the 
origin of Vela Jr. 

5. Conclusions 

We have presented our attempt to detect unshocked ejecta 
of the young SNR Vela Jr. by analyzing broad Ca II absorp- 
tion lines in spectra of background stars. We obtained and 
analyzed spectra of 14 stars across Vela Jr. using standard 
methods of spectral synthesis. We concluded that broad ab- 
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Fig. 8. Profile of 1.16 MeV line of ''''Ti. Panel a: profiles for 
SN Ibc (thin line) and SN Ic(h) (thick line) with spherically- 
symmetric distribution of ''^Ti. Panel b: profiles for SN Ic(h) 
with 44^j distributed in the external parts of bi-polar jets at 
inclination angles of 30° (solid line) and 60° (dashed line). 



sorption lines are absent. The 3cr upper limit to the depth 
of broad absorption lines is 0.04. We detected both low ve- 
locity and high velocity interstellar Ca II absorptions. The 
latter are attributed to the cold gas of radiative shocks 
propagating in clouds engulfed by the old Vela SNR. 

The absence of broad Ca II absorption lines and the 
constraints imposed by the flux of the 44^^^ gamma-ray line 
and the angular size of the SNR imply that only SN Ibc or 
energetic SN Ic (hypernovae) could have produced Vela Jr., 
if our estimates of the age and distance are correct. The ad- 
ditional constraint provided by the width of the 1.16 MeV 
'''^Ti line also supports the hypernova scenario for Vela Jr. 
origin. However, we emphasize the need for more reliable 
data on the 44 gamma-ray line profile and higher reso- 
lution imaging of Vela Jr. in the gamma line to verify the 
hypernova scenario. 
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